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The flapping flag instability occurs when a flexible cantilevered plate is immersed in a 
uniform airflow. To this day, the nonlinear aspects of this aeroelastic instability are largely 
unknown. In particular, experiments in the literature all report a large hysteresis loop, 
while the bifurcation in numerical simulations is either supercritical or subcritical with a 
small hysteresis loop. In this paper, this discrepancy is addressed. First weakly nonlinear 
stability analyses are conducted in the slender-body and two-dimensional limits, and 
second new experiments are performed with flat and curved plates. The discrepancy is 
attributed to inevitable planeity defects of the plates in the experiments. 



1. Introduction 



The scientific interest in the flapping flag instability dates back to Lord Rayleigh ( 1879 ) 



As a side remark in his famous paper on the stability of jets, he showed that an infinite 
membrane placed in an airflow is always unstable. Of course, the problem becomes more 
complex when bending rigidity and finite plate dimensions are taken into account (see 



Paidoussis: 2004 S helley fc Zhang||2011 for recent reviews) 

Theoretical models of this instability can be divided into two categories according to 
the plate aspect ratio. For large aspect ratios, a two-dimensional analysis is relevant, 



and Kornecki et al. ( 1976 ) were the first to show that the flow around the plate can 



be modelled using unsteady airfoil theory When the plate aspect ratio is asymptotically 
small however, the aerodynamic forces can be modelled using slender-body theory (iDatta 



& Gottenberg 1975 Lemaitre et al. 2005|. The studies in these two asymptotic limits have 
been recently generalised by Eloy et al. (2007) and Doare et al. (2011) who considered 



intermediate aspect ratios and confinement effects. 

The aforementioned stability analyses are all linear though and, to our knowledge, the 
nonlinear dynamics has never been addressed theoretically. This is one of the goal of this 
paper, the other being to compare these new weakly nonlinear analyses to experiments. 



Since the pioneering study of iTaneda (1968), different groups have performed exper- 



iments on the flag instability either for small aspect ratios (Datta & Gottenberg 1975 
Lemaitre et aL|2005 ) or for moderate to large aspect ratios ( Zhang et a/.|2000 Tang et al. 



2003[ |Eloy et aL|2008 among others). In this latter large hysteresis is always ob 

served at threshold or, said differently, the motionless state and the flapping state coexist 
and are both stable over a large range of airflow velocities. To quantify the importance 
of this hysteresis, we will refer to the hysteresis loop, (U c — Ud)/U c , where [Ud U c ] is 
the velocity range of bistability. In large aspect ratio experiments, this hysteresis loop is 
typically of 20%. For aspect ratios smaller than one however, hysteresis can disappear 
flEloy et a*.||2008l ). 
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Parameter 


Symbol 


Value 


Unit 


Plate length 


L 


3-30 


cm 


Plate height 


H 


1-20 


cm 


Plate surface density 


m 


0.25 


kgm" 2 


Plate bending rigidity 


D 


3.0 x 10" 3 


Nm 


Fluid damping 


V 


0.044 


kgm" 2 s 


Internal damping 


/' 


1.7 x 10~ 6 


kg m 2 s~ 


Wind velocity 


V 


0-65 


ms" 1 


Air density 


P 


1.20 


kgm" 3 



Table 1 . Characteristics of the experiments. 



For small aspect ratios, and negecting the nonlinearities of the aerodynamical forces, 



Yadykin et al. (2001) showed that elastic and inertial nonlinearities yield no hysteresis. 



In the two-dimensional limit, most numerical simulations have exhibited bistability, both 



for inviscid flows modelled with vortex methods (Alben & Shelley||2008 Michelin et al. 



20081 or for viscous flows modelled with Navier Stokes solvers (Zhu & Peskin 2002 



Connell k. Yue|2007 1 . The only exceptions are the studies of Tang et al. ( 2003 ) who did 



not consider aerodynamical nonlinearities and Tang & Pai'doussis ( 2007 1 who considered 
relatively short and heavy plates. Even when bistability is present, the hysteresis loop is 
much smaller than in the experiments: Alben & Shelley (2008) and Michelin et al. (2008) 
report an hysteresis loop of 2.5-4% and 4.5% while experiments in the same range of 
parameters exhibit loops of 20% or more. Several effects have been advanced to account 
for this discrepancy: blockage and confinement effects ( Tang fc Paidoussis|2007 ), planeity 
defects of the plate ( Tang fc Paidoussis|2007| |Eloy et a/.|2008 ), or damping effects (Alben 
fc Shelley||2008| . But, until now, none of these hypotheses have been tested 



This paper is organised as follows: in §[2] the experimental setup is briefly described; 

3 the physical model is introduced and the weakly nonlinear analyses are carried out; 

4 the experimental results are presented and compared to the theoretical predictions; 
and finally in §[5] these results are discussed. 



2. Experiments 

The experimental set-up is illustrated in figure [IJa. Experiments are performed in a 
horizontal low-turbulence wind tunnel of 80 x 80 cm 2 cross section. The flexible rectan- 
gular plates are cut from Mylar sheets whose physical characteristics are given in table [1] 
The bending rigidity D has been measured through deflection tests under gravity, and 
the fluid and material damping coefficients v and /i have been evaluated by measuring 
the damping of clamped plates in air at rest. 

The same protocol has been followed in all experiments and will be briefly described. 
Once the plate is clamped into the mast, the flow velocity is gradually increased. At 
small velocities the plate is motionless. Eventually, for a critical flow velocity U c , the 
plate starts to flutter with a large amplitude and a well-defined frequency. When the 
flow velocity is decreased, the plate returns to its stable state again at a different critical 
velocity Ud ^U c , thus leading to a hysteretic cycle. 

The motion of the plate is recorded with a CCD laser displacement sensor of spatial 
and frequencies accuracies of 1 /xm and 10 000 Hz. In the present study, the deflection was 
measured at the plate center. At this point and near threshold the deflection is always 
harmonic in time, and the amplitude A will refer to the peak amplitude of the deflection. 
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(a) streamlined mast 




Mylar plate 




laser sensor 



Figure 1. (a) Sketch of the experimental setup; (6) definition of the coordinates and plate 
dimensions; and (c) view of the one dimensional plate motion x(s,£) = (x,y). 



3. Modelling 

Consider a flexible plate of length L and span H clamped into a mast parallel to the 
vertical axis Oz and immersed in a uniform flow of velocity U (figure [lj)). Assuming the 
motion is independent of z, this plate obeys the Euler-Bernoulli beam equation 

mdfx + Dd*x - d s ((T)d s x) + vd t x + pd t dfx + (p)n = 0, (3.1) 

with s the curvilinear coordinate, x(s,t) = (x,y) the plate position, m its mass per 
unit area, D its bending rigidity, T the generalised tension in the plate that enforces 
inextensibility (i.e. ||<9 s x|| 2 = f), and the chevrons (.) denote the average along the span. 



The terms proportional to v and p in (3.f ) model the dissipation due to the fluid viscosity 
and the Kelvin- Voigt structural damping respectively (Pai'doussis 2004). The last term 
in (3.1) is due to the pressure jump across the plate, p(s, z,t), and n is the unit vector 
normal to the plate (figure [l]c) . In addition the deflection y(s,t) satisfies clamped-free 
boundary conditions: y = d s y — in s = 0, and d^y — d^y = T = in s = L. 

The Reynolds number being large, we will further assume that the flow is potential 
such that the total flow velocity is given by u = U + V0, with U the imposed airflow 
and <j) the perturbation potential. The potential <f> can be found by solving a Laplace 
problem with Neumann boundary conditions ensuring the impermeability of the plate 

Acj) = 0, with V4> ■ n = (<9 t x - U) • n on S, (3.2) 

where S is the moving plate surface. From <p the pressure jump can be calculated using 
the unsteady Bernoulli equation. 

Using L and L/U as characteristic length and time, the control parameters are reduced 
to three dimensionless numbers noted with asterisks: the reduced velocity U* , the mass 
ratio M* , and the aspect ratio H* 

U* = ULy/m/D, M* = pL/m, H* = H/L, 

and the dissipative properties are characterised by two coefficients 



vp m 2 D 2, p* — pp m ?D 



(3.3a-c) 
(3.4a, 6) 



where p is the fluid density. With this non-dimcnsionalisation the damping coefficients 
depend neither on the plate dimensions H and L nor on the flow velocity U and thus are 
constant in the present study : v* = 0.071 ± 0.001 and p* = 0.0014 ± 0.0003 (with 95% 
confidence intervals). 

Equations (3.1 3.2|) describe the nonlinear, fully coupled, fluid-structure interaction. 



We will now consider three approximations of these equations: the linear limit to address 
the instability threshold, and the weakly nonlinear limit both for an elongated plate and 
for an infinite plate, to study the nature of the bifurcation. 
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3.1. Linear model 

Assuming small lateral deflection (y <C L and d x y <C 1), the system (3.1-3.2) is linearised 
md 2 t y + Dd*y + vd t y + pd t d 4 x y + (p) = 0, 



A<^> = 0, with d y 4> = dty + Ud x y for y = and < a; < L. 



(3.5a) 
(3.56) 



This system is solved by the same method as in our previous papers (Eloy et al. 2007 



2008), except that dissipative terms are retained in the analysis. The main steps of this 



linear stability analysis are as follow. First, a complex angular frequency u> is assumed and 
the deflection y is expanded on Galerkin modes that satisfy the clamped-free boundary 
conditions. Second, the Laplace problem ( |3.5| ?>) for the perturbation potential is solved 
for each Galerkin mode in three dimensions, and the associated average pressure jump 
(p) is calculated. Finally, the eigenvalue problem obtained from the equation of motion 
( |3.5| a) is solved to obtain the global modes and their complex frequencies. If one of these 
complex frequencies has a negative imaginary part, the plate is unstable. 

This linear stability analysis allows to predict the critical velocity U*, above which the 
system is unstable, as a function of the dimensionless parameters M* , H* , u*, fx*. 

3.2. Slender-body nonlinear model 

In the slender body limit (i.e. H -C L), the aerodynamic force on the plate can be 
decomposed into two terms: a reactive force originating from the added mass of air 
accelerated when the plate moves, and a resistive force modelling the drag on the plate due 
to crossflow. As a result, the average pressure jump is decomposed as (p) = p rca c. +Presis., 
where p re ac. is the reactive part and has been calculated by Light hill| ( 197ij ) in the context 
of fish locomotion, and p les is. is the resistive part 



Pr. 



M (w — (uw)' 



Pr 



\pC d \w\w, 



(3.6a, b) 



where dots and primes denote differentiation with respect to t and s, u and w are the 
longitudinal and normal component of the plate velocity relative to the uniform airflow 
(such that x — U = ut + wh) 



u = —Ux + ; 



yy 



= Uy' 



xy + yx 



(3.7a, 6) 



Al = irpH/4 is the added mass of air, n = y" / (1 — y' 2 ) 1 ^ 2 is the plate curvature and Cd 



is a drag coefficient taken to be Cd — 1.8 for a plate (see |Buchak et aI1[20To| . 

Inserting (3.6-13.7) into (3.1 1 and projecting onto x and y give two coupled dynamical 



equations for x(s,t) and y(s,t). Following Yadykin et al. (2001 1, these equations are 
decoupled by first using the x-projection to eliminate the average tension (T). Then 
x(s,t) and its derivatives are eliminated by using the inextensibility condition. Finally 
the terms of order larger than y 3 are discarded yielding a weakly nonlinear dynamical 
equation for y(s,t) that can be expressed as 



C(y) + mf m (y) + Dfn(y) + Mf ieac .(y) + P C d f resis .{y) = 0, 
where £ is a linear differential operator on y 

C(y) = my + Dy"" + M (U 2 y" + 2Uy' + y) , 
and the other terms are 0(y 3 ) nonlinear terms 



f m (y) 
f D (y) 



y 



y""y' 2 



{y 12 + y'y') ds - y'> 

ly'y"y"' + y"\ 



y'y 1 



dsd 



s, 



(3.8) 
(3.9) 

(3.10a) 
(3.106) 



/reac. {y) 



y y 
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1 

y 



2(Uy" + H 
1 



U{y'y' 

S 

y'y'ds 



3y"y'y) - Zy'y'y 



II -2 

:V y 



y" / y' {u 2 y" + 2U y ' + y) ds, 



f IC s is .(y)^^\u y ' + y\(Uy' + y), 



y'y) ds + 
(3.10c) 
(3.10d) 



except /resis. which is of order y 2 . For the sake of brevity dissipative terms have been 
omitted in this analysis (we have checked that they were indeed negligible). 
The deflection y is now assumed to be of the form 

y{s,t) = (Ah {s) + /i 1 (s))e iw * + c.c. (3.11) 

where A is a small complex amplitude (\A\ <C L), ho is the solution of the linear problem 
£(/ioe Iwt ) = with clamped-free boundary conditions and the proper normalisation to 
allow comparison with experiments (i.e. h (L/2) = 1/2, such that y(L/2,t) = Acosuit), 
h\ gathers the terms of order greater than A, and 'c.c' stands for 'complex conjugate'. 

Using standard methods of perturbation theory, the weakly nonlinear amplitude equa- 
tion for A is found. It consists in finding and solving the linear adjoint problem, inserting 



the decomposition (3.11 ) into (3.8 1, and forming the scalar product with the adjoint solu- 



tion. It yields the following amplitude equation given in dimensionless form for A* = A/L 



w A 



ioj*U* Cl A 



+ (u* 2 c 2 + c 3 ) A* + (U* 2 c 4 + iu*U*c b )\A*\A* + 

'u* 2 CQ + iuj*U*c 7 -Lo* 2 c % +c 9 ) \A*\ 2 A* =0, (3.12) 



where w* = ujL 2 ^/mD, and the complex coefficients Cj depend on the control parameters 
M* and H*. Once the c,-'s are calculated, the solutions of the second degree equation 



(3.121 for uj* can be obtained. Then, values of \A* \ such that the imaginary part of uj* is 
zero give the saturated amplitude A*(U*). 



The first nonlinear effect comes from the resistive force (3.10 d) and gives the term 
proportional to |A*|^4* in (3.12). As a result, A* is proportional to (U* — U*) close to 
threshold (see figure^). This behaviour has to be contrasted with the usual pitchfork 
bifurcation for which A* is proportional to (U — U*) 1 ^ 2 . In any case, the bifurcation is 
not subcritical, and one therefore expect no hysteresis in the slender body limit. 

3.3. Two-dimensional nonlinear model 

For plates of large aspect ratios (H L), substantial assumptions have to be made 
to solve the weakly nonlinear problem analytically. First the plate span and length are 
assumed infinite. The problem then becomes two-dimensional, and it is further assumed 
that the deflection is a propagating wave of amplitude A 



y = A cos(kx — Lot) . 



(3.13) 



To calculate the flow around the plate, the perturbation potential is expanded in powers 
of A (i.e. <f> = A<j>i + A 2 cf>2 + A 3 (f>3 + •••), similarly to the methods used for the weakly 



nonlinear analysis of the Rayleigh-Taylor instability ( Nayfeh 1969 ) . Inserting this ex- 



pansion into (3.2) and using a Taylor expansion to evaluate <f> on the plate, the Laplace 



problem is solved for the first three orders yielding 

T(V -U) sin(fcx-wt) e Tky , 
= - \k{V - U) sin(2fcE - 2uot) e^ 2 ^, 



= T kk 2 (V-U) sin(kx-ut) e* ky + h.o.t. 



(3.14a) 
(3.146) 
(3.14c) 
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where V = tu/k is the wave speed of the deformation, the superscript ± corresponds to 
the upper and lower parts of the flow and 'h.o.t.' stands for 'higher order terms'. 

The pressure field is deduced from the potential 4> using the unsteady Bernoulli equa- 
tion. The y-component of the pressure force Fp is then calculated by evaluating the 
pressure jump using a Taylor expansion and projecting it onto the vertical axis. Keeping 
only the first harmonics and the terms up to order A 3 yields 



F P = F x (1 - | A 2 k 2 



), with Fi = 2p(V - U) 2 Akcos(kx - cut). 



(3.15) 



The nonlinearities originating from the inertial and elastic terms can be evaluated 



by using the expressions (3.10 a, b). Note that derivatives with respect to s have to be 
transformed into derivatives with respect to x using the chain rule. Note also that, when 



calculating (3.10 a), boundary terms arising from the integrations are neglected because 
y(x,t), as given by (3.13), do not satisfy the clampcd-free boundary condition. 

Neglecting the dissipative terms, the different nonlinearities can now be gathered to 
obtain a weakly nonlinear dispersion relation in dimcnsionless form (with A* = A/L and 

k* = kV) 



1 



1 



A* 2 k* 2 



1 



1 



A* 2 k* 2 



M* 



u*ky i 



lA* 2 k* 2 



o, 



8" " J " \ 2 J \k* 

'(3.16) 

where the first term corresponds to inertia, the second to the bending force, and the last 
to the pressure force. Solving this second degree equation for lo* and finding the values 
of A* for which the imaginary part of lj* is zero gives the saturated amplitude A*(U*). 

When the mass ratio and the wavenumber are the same as in the experiments (M* = 
0.6, k* Ri 3ir/2), the weakly nonlinear dispersion relation (3.16) predicts a subcritical 
bifurcation, with a hysteresis loop of approximately 0.1% (see figure^). The bifurcation 



is always subcritical, except for small values of the mass ratio, M* < fc*/10, for which 
it is supercritical. This prediction should be considered with caution because several 



approximations have been made to obtain the dispersion relation (3.16): in particular, 



the plate has been assumed infinite to calculate the pressure forces and we know from 
the linear stability analysis that finite size effects can be important. 

Yet, a comparison with the inviscid numerical simulations found in the literature give 
good qualitative agreement. For M* = 0.2, k* = 3tt/2, the bifurcation is found to 
be supercritical (since M* < fc*/10) as in | Tang fc Paidoussis| |2007[ ); for M* = 3.3, 
k* = 5tt/2, the hysteresis loop is equal to 5.2% which compares well with the 3.4% 
reported by|Alben & Shelley| (|2008|); finally for M* = 10, k* = 7tt/2, the hysteresis loop 



is found to be 12.8% while Michelin et al. (20081 found 4.5% 



4. Results 

To examine the dynamics near threshold, two sets of experiments have been carried 
out, either with flat plates, or curved plates (see the inset of figure (2)3) ■ Flat plates were 
cut from unused Mylar sheets and clamped into the mast, whereas curved plates were 
first heated on a curved surface for 10 s. In this latter case, an intrinsic curvature of the 
order of 0.1/1/ in the vertical direction subsisted (this curvature is exaggerated in the 
inset of figure [2)3). For both flat and curved plates, each experiments was reproduced five 
to ten times to assess the repeatability and to extract statistics. 

Experimental values of the flutter amplitude A*, measured for M* = 0.6 and three 
different values of H* , are plotted in figure [2] It first shows that the deflection A*{U*) 
does not depend on the intrinsic curvature of the plate and that the repeatability on 
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Figure 2. Flutter amplitude as a function of the velocity U* for a mass ratio M* — 0.6 and 
aspect ratios: H* = 0.5 (a), H* = 1.5 (b), and H* = 1.0 (c). Open and filled symbols correspond 
to experiments with flat plates and curved plates (see the inset in (a)). The corresponding 
average thresholds U2 and U* are represented with the same convention (exact values are given 
in table [2b. The lines in (a-c) represent the predictions of the weakly nonli near analys es in the 
slender-body (-) and two-dimensional limits (- -), as calculated from ( |3. 12 1 and ( |3.16[ ). Finally 
the square of the amplitude is plotted in (d) for three experiments with different aspect ratios 
as labelled (lines are best linear fits near threshold). 



this measurement is excellent. These plots also show that the weakly nonlinear analyses 
developed in §|3]give good predictions of the saturated flutter amplitude. In the elongated 
plate limit (figure [2)3), the analysis predicts a supercritical bifurcation with A* propor- 
tional to (U* — U*) near threshold. However, as illustrated in figure [2]i, A* can always 
be reasonably fitted as A* ~ (U* — U*) 1 / 2 for all aspect ratios, which indicates that the 



resistive force modelled by (3.10 d) may not be valid for small flutter amplitudes. Since 



this drag force is caused by vorticity shed from the top and bottom edges of the plate, it 
probably means that this shedding occurs only when the flutter amplitude is sufficiently 
large. 

Experimental data for H* = 1.5 are plotted in figure^ together with results of the 
two-dimensional analysis. In the inset of the figure, the scale of U* has been expanded to 
emphasize the subcritical nature of the bifurcation. This subcriticality does not appear 
on experimental data (see figure^), but this could be due to the moderate value of the 
aspect ratio. For H* = 1.0 (figure [2]c), the measured flutter amplitudes lie between the 
predictions of the slender-body and two-dimensional approximations. 

Statistics on the the critical velocities have also been performed and the results are 
reproduced in table[2](see also the arrows in figure [2] indicating the mean values of U* and 
Uj ). These statistics show that U% exhibits a good repeatability and has indistinguishable 
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H* = 0.5 H* = 1.0 H* = 1.5 

m u: u* d u* u* d u* c 

Flat plates 11.7 ±0.2 12.3 ± 0.3 9.4 ±0.1 12.0 ± 1.8 9.0 ±0.4 15.4 ± 4.9 

Curved plates 11.5 ±0.1 19.9 ±3.7 9.5 ±0.1 20.8 ± 2.9 8.7 ±0.4 12.8 ±1.5 

Table 2. Statistics of the measured critical velocities U£ and U* for M* = 0.6 and three 
aspect ratios H* = 0.5, 1.0, and 1.5. For each case, the first figure is the mean and the second 
the standard deviation measured on 5 to 10 similar experiments. 



values whether the plate is curved or flat. On the contrary, the mean value of U* changed 
drastically between flat and curved plates. In addition, the standard deviation for U* is 
consistently larger than for U*[, showing a relatively poor repeatability. 

These results show that U* is extremely sensitive to curvature. As it was shown by 
Peake (1997), the plate curvature can be taken into account in the linear regime by 



adding a spring foundation term to the equation of motion (3.1), with spring rate per 
unit surface, Ehn 2 , with E = 1.5 GPa the plate Young's modulus, h = 280 /im the plate 
thickness, and k ss O.lrn -1 the curvature. In dimensionless terms, this corresponds to 
adding a constant term of order 5 2 /h 2 to the dispersion relation (3.16), where S is the 
typical deflection due to curvature. For the curved plates used here, this dimensionless 
spring rate is approximately 200 and leads to a 20% increase of the instability thresh- 
old, in qualitative agreement with the experimental observations (this term is only an 
approximation because nonlinearities should be considered when S > 0(h), Audoly & 



Pomcau 2010). As soon as planeity defects are of the order of the plate thickness, and 



this is almost impossible to avoid in practice, one then expects a stiffening effect and a 
delayed instability. 

When the plate starts to flutter however, one expects this curvature to be 'ironed out', 
as originally hypothesised by Tang & Pai'doussis (20071. This ironing out is due to the 
prohibitive cost of having a non-zero Gauss curvature that would introduce stretching 
energy in addition to bending energy. 

In brief, we propose the following scenario. Because of curvature or planeity defects, the 
plate is stiffen and the instability is delayed. Once the plate flutters however, these defects 
are ironed out and the decreasing threshold U% is the one predicted by the analyses. This 
scenario is coherent with all the experimental observations. It explains why A* and U£ are 
identical when curvature is added, why U* has poor repeatability and can be increased 
by curvature. It also explains why the hysteresis loop tends to increase with the aspect 
ratio for flat plates (curvature defects should also have a zero Gauss curvature and thus 
are more likely to be in the vertical direction when H* ^ 1). Note that for the largest 
aspect ratio, H* = 1.5, the curved plates have a lower mean threshold U* than flat plates, 
but this is not statistically significant because of insufficient data. 

The present stiffening scenario might not be the only cause of hysteresis: in particular, 
it cannot explain the bistability observed in the soap film experiments of |Zhang et al. 



(2000). However, the alternative explanations found in the literature are ruled out by the 
present set of experiments: both blockage and damping effects would not lead to poor 
repeatability of U* and a drastic increase of U* when curvature is introduced. 

Finally, it should be noted that no amplitude smaller than A* w 0.06 could be mea- 
sured, for all aspect ratios and despite the care taken. It is thus difficult to give a definite 
answer on the nature of the bifurcation. Nevertheless, the amplitude clearly behaves as 
A* ~ (U* — U*) 1 / 2 near threshold (figure ^p) and we have shown that the hysteresis 
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Figure 3. Instability threshold as a function of H* for fixed M* = 0.6 (a) and as a function 
of M* for fixed H* = 0.5 (6) for flat plates. Open and filled triangles correspond to and 
respectively. The thick line is the result of the linear stability analysis detailed in §3.1| and, the 
thin line is the same analysis without dissipative terms (i.e. for v = \i = 0). The filled circles in 
(b) mark the departure from periodic motion and the dotted line indicates the location where 
the second and third eigenmodes have the same frequency. 



is mainly caused by planeity defects. Without these defects, we would thus expect a 
supercritical bifurcation, in agreement with the theory (the bifurcation could become 
subcritical for larger H* or larger M* , but then, plates may sag under their own weight). 

Figure [3] shows how the thresholds U* d and U* vary with aspect ratio for constant mass 
ratio, or vary with mass ratio for constant aspect ratio. This measurement is compared 
with the predictions of the linear stability analysis and the agreement is excellent for Uj , 
confirming the proposed scenario. It also shows that taking into account the internal and 
viscous damping gives a better prediction of the instability threshold. 

Figure[3^ also shows that, for sufficiently large airflow velocities, the flutter is no longer 
periodic. This departure from periodic motion occurs apparently when the unstable mode 
(with wavenumber k* ss 3n/2) has the same frequency as the next eigenmode (with 
k* 57r/2). A nonlinear interaction between these modes is thus likely to be the cause 
of this secondary bifurcation, which is currently being investigated. 



5. Conclusion 

In this paper, we addressed theoretically and experimentally the nonlinear dynamics 
of the flag instability near its threshold. Weakly nonlinear analyses have been carried out 
both in the slender body limit and in the two-dimensional limit. In the former case, the 
bifurcation is always supercritical, while it can be subcritical in the latter case for large 
enough mass ratios (in agreement with numerical simulations of Tang fc Pa'idoussis|2007 
Alben fc Shelley||2008[ |Michelin et aZ.||2008 l 



In the experiments, the dynamics near threshold is complicated by inherent planeity 
defects, which usually lead to large hysteresis. These defects are likely to be the main 
cause of discrepancies between numerical simulations and experiments. They are also 
responsible for the poor repeatability of the measured instability threshold when airflow 
velocity is increased. However, the threshold measured for decreasing velocity is both 
repeatable and in excellent agreement with the predictions of the linear stability analysis. 
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